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Recent results from Belle and BaBar Collaborations hint at a small sin2<^i, while the measured 
B — + TTTT rate also seems to be on the low side. Supersymmetric (SUSY) models with down squark 
mixings can account for the deficits in both cases. By studying the origin of SUSY contributions 
that could impact on B'^-B^ mixing and B — > tttt decay, we find that the former would most likely 
arise from left-left or right-right squark mixings, while the latter would come from left-right squark 
mixings. These two processes in general are not much correlated in the Minimum Supersymmetric 
Standard Model. If the smallness of _B ^ tttt is due to SUSY models, one would likely have large 
B py from chiral enhancement, and the rate could be within present experimental reach. Even 
if B ^ P7 is not greatly enhanced, it could have large mixing dependent CP violation. 



PACS numbers; 12.60.Jv, ll.30.Er, 13.25.Hw 
I. INTRODUCTION 

The CP asymmetries in B" J/ipKs, J/^Kl de- 
cays have been studied by several experimental groups 
|^|-^. It is well known that one of the phase angles of 
the Standard Model (SM) unitarity triangle, sin20i, can 
be measured via the asymmetry, 

_ r(B^>it) ^ j/i^Kl) - r(^o(t) ^ j/tpK"s) 

= — sin 201 sin Am i . ( 1 ) 

The CDF Collaboration finds sin2(/>i = 0.79 [§ with 
Tevatron Run-I data, while the OPAL and ALEPH Col- 
laborations give sin 201 = 3.2+^-^ ± 0.5 @, 0.84+°'^^ ± 
0.16 1^, respectively. Recently, however, the BaBar and 
Belle Collaborations announced their results on the mea- 
surement of this asymmetry. The Belle Collaboration 
reports sin20i = 0.58 +°-34+° °o g, while the BaBar 
Collaboration gives the even smaller sin20i — 0.34 ± 
0.20 ±0.05 When combined with previous CDF and 
LEP results, the average value is sin20i = 0.48±0.16. 
While this is consistent with the Cabibbo-Kobayashi- 
Maskawa (CKM) fit value of sin20i = 0.698±0.066 § or 
0.47 < sin 201 < 0.93 (95%C.L.) 0, the central value is 
rather small. It could be hinting at the presence of new 
physics effects, especially if the value persists. In this 
case, we may need a large new physics contribution ||] 
in B^-B^ mixing, comparable to the SM amplitude, to 
account for the smallness of aj/^Ks- This is because it is 
very hard for new physics to affect the Cabibbo favored 
h CCS decay amplitude. 

The first result on the charmless decay mode B^ 
TT+TT^ was given by the CLEO Collaboration, giving 
Br(B*' TT+TT-) = (4.3+}-^ ± 0.5) x lO^^ g. BaBar 
and Belle Collaborations also reported recently their re- 
sults, Br(BO TT+TT") = (4.1±1.0±0.7) x lO'^ |o|, 
(5.9+2;f ± 0.5) X 10-6 [|ii|, respectively. Note that 
the BaBar and Belle measurements are all lower than 
their reported results at summer 2000 conferences p2| , p^ . 
The combined result with averaged Br(i3'^ tt+tt^ j = 



4.4±0.9 seem to be on the low side when compared to the 
SM prediction using factorization approach, Br(i?'' 
TT+TT") ~ 10 X 10-6 ||l||, for 03 - 60°, and remains true 
when compared to the QCD factorization result |l^ of 
Br(S" TT+TT-) - 8 X 10-6. In a recent work on QCD 
factorization, Br(TT+TT-) can be close to the experimen- 
tal value, however, the lowness of the rate would subject 
the SM to considerable stress ||l6|. A 20%-40% or more 
reduction in the branching ratio is welcome. In SM, the 
tree amplitude dominates over the penguin amplitude, 
which is about 30% of the former. Thus we may need 
large contribution from New Physics if it is responsible 
for the smallness of the rate. 

We have two cases where we are in a situation that 
New Physics contributions should be large, if it is re- 
sponsible for the smallness of measurements. As one 
of the leading candidates for new physics, supersymme- 
try (SUSY) helps resolve many of the potential prob- 
lems that emerge when one extends beyond the SM, 
for example the gauge hierarchy problem, unification of 
SU(3)xSU(2)xU(l) gauge couplings, and so on In 
the context of SUSY, we then ask the following ques- 
tions : Is it possible for SUSY models to account for the 
smallness in both processes? If so, are they correlated, 
since both of them are 6 — > d flavor changing processes? 
Since New Physics contributions would be large, can we 
find other related effects? 

To analyse SUSY contributions, we follow the ap- 
proach of Ref. ||l8|. As will be discussed later, gluino ex- 
change diagrams induced by d-h mixings give dominant 
contributions in both of the above mention processes. We 
do not aim at constructing any explicit models, hence we 
have not considered other flavor changing processes such 
as K^-K^ mixing, D^-D^ mixing, Br(i3 Xg^) and 
the neutron electric dipole moment, etc., since these are 
controlled by other parameters. Our strategy has been 
simply to study the implications on d-b squark mixings 
from new data on Bd mixing and B tttt decay, and 
make inference on other modes, such as -B — + p7 which 
is quite correlated with effects in i? — > tttt. We have as- 
sumed that models can be constructed such that SUSY 
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can impact on the modes considered here, but do not run 
into trouble with other stringent low energy constraints 
(see e.g. Ref. for Bd mixing case). 

We organize this paper in the following way: We dis- 
cuss SUSY contribution to B^-B^ mixing, B ^ tttt and 
radiative B decays in the next two sections. We then 
give some discussion, followed by conclusion in the last 
section. 



II. B^-B'^ MIXING IN SUSY MODELS 

The effective Hamiltonian for B^~B^ mixing from 
SUSY contributions is given by 

-Y,C,0,, (2) 



where, 
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together with three other operators 01.2,3 (and associ- 
ated coefficients Ci) that are chiral conjugates (L ^ R) of 
C'i,2,3- There are contributions from gluino, neutralino, 
charged Higgs and chargino exchange diagrams jl^ . We 
note that, due to Majorana property of the gluino, gluino 
exchange diagrams can be divided into the usual box di- 
agram and the so called crossed diagram. By using the 
double line notation of 't Hooft , it is easy to see that 
the former has a color Nc factor, while the latter does 
not. In general, therefore, the leading SUSY contribu- 
tion comes from gluino box diagrams, where we have 
and Nc enhancement, although it is possible that in some 
parameter space, such as small tan/3 and when superpar- 
ticles are light, charged Higgs and chargino contributions 
may become important |2l| . 

It is customary to take squarks as almost degenerate 
at scale to. In the following, we give contributions from 
gluino exchange diagrams and make use of the mass in- 
sertion approximation 22|. In quark mass basis, one 
defines 0, 



"qAB 



('^g)!4s/"^^ 



(4) 



which is roughly the squark mixing angle, are squark 
mass matrices. A, B ^ L,R, and i,j are generation in- 
dices. For notational simplicity, we shall suppress in what 
follows the index pair ij (13 for a d-b mixing angle) as 
well as the subscript d. 

The gluino exchange contributions to the Wilson coef- 
ficients are [|l8|. 
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where Xgq = m^/rn^, and (71^2.3 are given by changing 
L i? in Ci_2,3, respectively. As noted before, the terms 
containing N^ are from the box diagrams, while those 
containing 1 are from crossed diagrams and l/-^c 
are from subleading terms of box and crossed diagrams. 
The loop functions are given by 



'f&{^gq)i ^gqfQ{Xgq)\ 



dm'^dm^ 



dk' 



(fc2 -I- TO?)2 

[e, ml] 



(6) 



(fc2 + m?){k'^ + to'2) 



m' — 'Tri 



which agree with Ref. ||l8[. Note that feixgq) is always 
positive, while feixgq) is always negetive. It is useful to 
give the asymptotic forms of these functions, 

[ ~ f&{Xgq)^ XgqfQ(Xgq)\ 

l/(3a;|q), for m < (xgg > 1), 

1/30, 1/20], for m~g ^ to {xgq = 1), (7) 

1/3, -Xgq\nxgq\, for rUg <t:m [xgq < 1). 

By using these asymptotic forms, it is easy to show that, 
\fti{xgq)\ < \xgqf6{xgq)\ for fh < TOg and vice versa. 
Therefore, one must have a zero in Ci (Ci) for some 



Parameter 


Value 


Parameter 


Value 


P 


0.169 


V 


0.362 




230 ± 40 MeV 


mB 


5.2794 GeV 


TjB 


0.55 


mt(mt) 


170 GeV 


fJ-B 


2.5 GeV 


mbifiB) 


4.88 GeV 




0.276 


MSUSY 





TABLE I. The input parameters used in this section. 



2 



value of Xgq hence a sign change when passing through 
it. One can also show that the cancellation is between 
the above mentioned two class of diagrams. 

After obtaining these Wilson coefficients at SUSY scale 
MgusY, we apply renormalization group running to ob- 
tain B mass scale values. The renormalization group 
running of these Wilson coefficients including leading or- 
der QCD corrections is given in Ref. Since B^-B^ 
mixing is a AB — 2 process, we need a power of S in 
each of the internal squark lines to change flavor. There 
are altogether six combinations 

and SlrSrl, as is evident from Eq. (I 

To obtain Am^^, we use Attlb^ — 'A^ri 



Bd mixings 
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and 
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where is the SM contribution, its value is well 

known [EJ] . The vacuum insertion matrix elements of Oi 
are given in Ref. p^ . These matrix elements are mod- 
ified by bag-factors to include non-factorizable effects. 
For simplicity, we assume the bag-factors for matrix el- 
ements of 02-5 are equal to Bb^, which is calculated 
for Ol. In the subsequent numerical analysis, we take 
fs^Bs^ = (230 ± 40) MeV For CKM matrix ele- 

ments, we take |X^b/AVcfc| — 0.41 and ^3 = 65°, hence 
\Vtd\ X 10^ = 8.0 to get Am|^ ~ 0.54 ps-\ which is close 
to the experimental value of Am^^ = 0.484 ± 0.010 ps~^ 
psf . We summarize input parameters used in Table I. 

In F igs. and ||, we show estimated limits of all six 
y^\6 S\s over the parameter space of sub-TeV gluino and 
squarks. The limits are taken such that the SUSY 
contribution to Bj mixing matrix element, IM12I, is com- 
parable with the SM result, hence large interference ef- 
fects could in principle occur that can give low aj/^^^. 

Note that by assuming the same bag-factor for Oi, the 

^1/2 

uncertainty on fs^B^^ does not show up in these fig- 
ures. Squark mixing angles that are much larger than 
those shown would give too large a contribution to Bd 
mixing and would require fine tuning to satisfy the exper- 
imental result. For mixing angles that are much smaller 
than those shown, they will not be able to generate 
large enough interference effect to reduce the asymmetry. 
Therefore, the limits shown in these figures may serve as 
upper limits from Am^^ constraint on one hand, and 
serve as roughly the required values to give impact on 

From Fi gs. B] and S w e see tha t the limits on 
\/\SllSrr\, \Slr,rl\ and ^/\SLRSRL\ are all of order 
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FIG. 1. Limits on 5ll,rr and SllSrr obtained by assum- 
ing IMfa^s^l < IMfa^^l in mixing. 
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FIG. 2. Limits on Slr.rl and SlrSrl obtained by assum- 
ing IMfaUS^I < |Mff I in mixing. 



few %, with \/\6llSrr\ as the most sensitive source 
for Bd mixing. This can be understood from Eq. (||), 
where \/\SllSrr\ in C4 has the largest Nc factor, while 
there is also RG enhancement Furthermore, we 

see from Eq. (||) that the dominant Sll^rr and 
terms are proportional to Xgqfe{xgq)/rh?, while domi- 
nant (5^,fl(5i^L terms is proportional to /e(a;gq)/m^. There- 
fore, the bounds on ^/\SLLSRR\ and \Slr,rl\ are roughly 
proportional to y/rh'^ /[xgqfQ{xgq)], while the bound on 

\/\6lrSrl\ is roughly proportional to yjm'^ /\fe{xgq)\, 
such that Figs. |l|(b) and||(a) show similar behavior that 
is different from Fig. |^(b). 

The order of magnitude of these figures can be un- 
derstood by a simple dimensional analysis. Compar- 
ing the SUSY vs SM box diagram contributions, we 
find S^< VtbV*a{aw/as){mtm/M^)il/VN,) < 2.8 x 
10^^ (to/500 GeV), which is very close to the limits on 
Slr,rl\, \/\SlrSrl\ and ^/\hJARR\ from Figs. |^(b) and 
However, the limit on as shown in Fig. |^(a) 

are of order few 10% and do not obey this estimation. 
This rather different behavior is because of the possible 
cancellation between Xgqfe{xgq) and feixgq) in Ci(Ci), 
which can weaken the bounds. A total cancellation is 
reflected in the valley along Xg^q ~ 2.43 where \Sll,rr\ 
is not constrained by \Mf^^^\ ~ |M™|. From Eqs.' (|) 
and (^, we see that Ci is dominated by Xgqfe{xgq) for 
Xgq 3> 1 and dominated by fe{xgq) for Xgq <C 1. There- 
fore, appart from the distortion due to the cancellation 
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discussed earlier, the upper left part of Fig. ||(a) is simi- 
lar to Fig. |l|(b) and Fig. ||(a), while the lower right part 
is similar to Fig. ||(b). Note that whenever we obtain 
a bound that is greater than 0(1) in the squark mixing 
angle (5, it should be interpreted as signaling the need of 
a large squark mass splitting which invalidates the ap- 
proximation of Eq. (||). 

It is clear that, to obtain low clj/^Ks^ need suit- 
able SUSY phase to have destructive interference with 
SM. But the minimum requirement is the SUSY ampli- 
tude should be large enough to allow for such an inter- 
ference effect. As shown in this section, it is possible for 
SUSY models to give Bd mixing that is comparable to 
SM result hence lead to large interference effect. Mixing 
angles of few % in left-right squark mixings or few % to 
few 10% in left-left, right-right squark mixings are suf- 
ficient to achieve this. The case with both left-left and 
right-right mixings {5ll5rr) is most sensitive to mixing 
angles, while left-left or right-right mixing alone are the 
least sensitive, and could even be totally insensitive for a 
fine-tuned parameter space near Xgq = ra-^/rh? ^ 2.43. 

For illustration, we pick a point from Figs. ||, ||, say 
(m, rUg) = (500, 400) GeV, and show the SUSY effects 
on Bd mixing. For the LL-RR mixing case, the bound is 
0.013, as shown in Fig. |l|(b). We show in Fig. 0, Ams^ 
and sin2$B_^, induced by 5rr = \5ll\ = (30%, 50%, 
70%) X 0.013, vs. aig{5LL), respectively. The a/| 5ll5rr\ 
are (30%, 50%, 70%) of the bound shown in Fig. |l](b) 
for the particular m, nig. Larger oscillating amplitudes 
in the figures correspond to larger S^l^rr. Similar re- 
sults will be obtained by using (30%, 50%, 70%) of the 
corresponding bounds of other points on (m, rrig) plane, 
shown in Fig. ^ and Fig. ^. The horizontal band in 
the left (right) figure is 2ct(1(t) experimental range. The 
uncertainty of the predicted Arris^ is due to the ^ 17% 

uncertainty of fsd^d shown in Table I. This fac- 
tor does not enter arg(M]^), and thus sin2$B^ shown in 
Fig. H (b), as we assume all bag- factors for different Oi 
are the same. By using a (5 at 30% of the bound value, 
Bd mixing does not differ much from the SM prediction, 
while for the 50% case, it starts to show interesting de- 
viation with sin2$B<j as low as 0.53. Note that the SM 
gives sin 201 = 0.73 by using our input parameters. For 
a larger S, such as 70% of the bound, it can further lower 



s'm2^Bd to 0.3. Although, in this case not all B,rg{SLL) 
are allowed, due to the Ato^^ constraint, we still have 
plenty of allowed region for this phase. From these fig- 
ures, we see that by using 50%-70% of the 5 bounds, low 
sin 2<i> can be easily obtained. 



III. B ^-KTT AND p7 DECAYS IN SUSY MODELS 

The effective Hamiltonian for charmless b ^ d decays 
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KbKd(ciOi + C2O2) - VM c,o.. 
-vMiCgdg + c'gd'g) 



4 = 3 



(10) 



where, as a matter of convention, we factor out a CKM 
factor VtbV^^ even for SUSY contributions. The operators 
are defined as, 



01 = 

02 = 

03(5) = 
04(6) = 

07(9) = 
08(10) = 



ujfj^Lbd'y'^Lu, 

d-f^Lbqj''L{R)q, 
da-fi_,Lbi3 qpj^'L{R)qa, 

^dj^LbQgq^RWq, 
3 - 

-dal^Lbp Qqqf)YR{L)qa, 



(11) 
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^Ma,,T'^'^^R{L)b q^T'^q, 
An q^ 



where — (1^ 75)72, Og"* arises from the dimen- 

sion 5 color dipole operator, and q = pb ^ Pd- Note that 
with New Physics, one may also have chiral conjugates of 
Oi_io with Wilson coefficients defined as c[_iq. The Wil- 
son coefficient Cg'^ and color dipole operators are defined 
in the effective Hamiltonian for b — > ^7, dg transitions, 

H^s = - ^"^'^f^ ie d [C, R + C'^L] a.^F^'^b 

+g d [Cg R + C'g L] a,,,T"Gr&} , (12) 



Parameter 


Value 


Parameter 


Value 


Fo(0) 


0.30 ± 0.04 




133 MeV 


n-iTT 


140 MeV 


TB 


1.548 ps 


niuifJ-B) 


2 MeV 




4 MeV 


const. 


0.2 GeV 




0.5 GeV 




1.5 GeV 




m^/3 



TABLE II. The input parameters used in this section. 
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ai 


d 


6 d 


ai 


b ^ d 


b ^ d 


ai 


1.0463 


1.0463 


a.0 


-0.0473 + 0.009H 


-0.0674 + 0.0029i 


0.2 


0.0502 


0.0502 


0.7 


0.0002 


-0.0000 


a-i 


0.0049 


0.0049 


OS 


0.0005 


0.0004 


a4 


-0.0337 + 0.0091i 


-0.0538 + 0.0029i 


acj 


-0.0093 


-0.0095 


a.5 


-0.0044 


-0.0044 


aia 


-0.0013 


-0.0014 



TABLE III. The OiS in SM for b ^ d and 6 ^ d for tjjs = 65°. 



where we have neglected irid, C^^g — + C^^f^ are 
the sum of SM and new physics contributions, while C'^ ^ 
come purely from New Physics. 

For B'^ — > TT+TT^, using factorization approach, we find 
01111 
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where Aa, is defined as a, 



with 
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(15) 



i/Nc, 



for even (odd) i. Using input parameters shown in Ta- 
ble II, the chiral factor R = 1.33. The expression for 
Stttt is somewhat different from the one given in Ref. ||2^] 
because of the treatment of q'^ in •* . We have used 

diaf,^ q" T^Rh qj^'T'^q 

= -(J7^ ^bT^Rb + d ^dli^T-Rb) {qi^T'^q) 
+dT''Rhpb^,[q{2 ^b- 4)rTy]. (16) 



where the ^ term can be dropped because of current 
conservation. By using heavy quark symmetry, pb —^pb, 
it is then straightforward to use factorization approach 
to obtain St^tt 



given in Eq. (15) 



For the form fac- 
tors involved, we have the relations /"""(to^) — 
F.iml), riml) ^ [ml/ml - l)[Fo(^^) - F^{ml)\ 
and 4mbh{ml) = ]+{ml) - /"(m^). Using Fo(0) = 
i^i(O) — 0.30 ± 0.04 and monopole form factors for Fo.i, 
with pole masses given in Ref. iQ, it is easy to show 



that /+(m2)/Fo(m2) 
4mb h{m.l)/ Foiml) - 



^ 1, r{ml)/Foiml) ~ 0, and 
1. Therefore, 



1 fn 



2N? 



-0.76, 



(17) 



which is not far from the value of —0.80 computed from 
Eq. (^5|), and also close to the value given in Ref. [p8| . 
Note that it is insensitive to Nc and the chiral factor R. 
The opposite sign between (a^, Cg) and (a-, C'g) can be 
easily understood by using parity transformation. We 
note that the color dipole term is sensitive to {q'^), which 
is usually taken to be between mf/i and ml/2. We use 
(g^) Tnl/i- The dependence will affect the color 
dipole contribution by ±33% at amplitude level. 

Using input parameters shown in Tables ||, ^ and 
following the approach in Ref. [|4j, we obtained the nu- 
merical values of a^s in SM, as shown in Table In 
SM, Cf^ = -0.15 while Cf^ is highly suppressed by 
the V— A nature of weak interaction and the smallness 
of md- The contribution is about 3% of tree ampli- 
tude. Under factorization approximation, the SM gives 
Br(i3° TT+TT-) - 10 X 10"^ with - 3.3% asymmetry, 
as defined as, 



Q.-rr+'Tr 



Br(BO ^ TT+TT") - Br(BO tt+tt-) 
Br(BO TT+TT-) -I- Br(BO tt+tt") ' 



(18) 



Since this process is tree dominant, we need large con- 
tribution if the rate is to be reduced by New Physics 
contribution. 

In SUSY models, we can have gluino, neutralino, 
chargino and charged Higgs exchange contributions to 
B — > TTTT. Because of Nc enhancement and different sen- 
sitivities of photonic vs gluonic penguins ||30[| , and since 
one does not suffer from Br(i3 Xs^) constraint, we 
see that gluino exchange gives dominant and interesting 
contribution to _B — s- tttt compared to other superparti- 
cles. There are two types of diagrams: The gluino box 
and the gluino penguin. The former as well as the Fi 
term (the quark chirality conserving vertex term) of the 



*When compare to the QCD factorization approach, we may 
use effectively (g^) ~ rni/Z - m\/2A in the color dipole con- 
tribution |29|. 
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latter contribute to a^, and only depend on one power 
of Sll.rr- The F2 term (the quark chirality flipped ver- 
tex term) of the gluino penguin contributes through Cg''' 
with all types of squark mixings. We use similar formulas 
as in Ref. ||T^ for SUSY contribution and those of Rcf. 
for RG running. 
The gluino box and Fi gluino penguin give, 
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(19) 
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where c[ are obtained from Ci by replacing L R, Bi{x) 
are from gluino box, P{x) = [C2(G)/2 - C2(i?)] Pi (a;) + 
G2{G)P2{x)/2 are from the Pi-term, and C2{G) = 
and G2{R) = [N^ - l)/2Nc are Casimirs. The Pi(2){xm) 
term corresponds to gluon attached to squark (gluino) 
line. The gluino box contribution to 03^4 (5^6) ^-re due 
to di-fei mixings in one of the squark line, while the 
other squark line is qL{R}- The leading terms in Nc are 
the Bi{x) and B2{x) terms of C4 and cg, respectively. 
Explicit forms of Bi{x), Pi{x) can be found in Ref. [[l8[ . 
They can be expressed as, 



[-4Bl{Xgg), B2{Xgg)] 

d 



9m'2 



(fc2 + m?)2(fc2 + m'2) 



(20) 



[k 



2 m?] 



(fc2 + m?) 



&[Pl{Xgq), -P2{Xm)\ 

d 



dm? 



dk' 



(fc2 + mj^){k^ + m?) 

2fc2- 



g' 



3ml 



(fc2 + m2)3 ' (fc2 + r71?)3 



(21) 



Note that Bi{x), P2{x) are always positive, while B2{x), 
Pi{x) are always negative. It is useful to give the asymp- 
totic form of these functions, 



[Bi(x), ^2(2;), P{x)] 

"lna;/j;2, -2/a;, -l/(9a;)], x > 1, 

1/12, -1/3, 8/45], X = 1, (22) 

1/2, 4xlnx, -21na;] , a; < 1. 



G. 



The gluonic and photonic penguins are closely related 
as for C^'^ and C^'^ 

nas Qd2C2{R) 



The formulas for Cj^ and Cg'^ from gluino exchange are, 



V2GFVtbV, 



td 



SLLg2iXgq) 

SLRg4{Xgq)}, (23) 



mb 



(Jll \2G2{R)-G2(.G)\g2{xgq) 



-G2{G)g^{x~gq) 



+ ^ ^Lfl [C2(G) - 2G2(i?)] g4{xgq) 
mb 



+C2(G)g3(xg,-)j|, (24) 

with the chirality partners G^ g obtainable by interchang- 
ing L i? in the 5's, is the electric charge of the down 
type quarks and the functions gi{xgq) = —m^ (d/dm^) 
[Fi{xgq)/m'^], where Fi{x) are given in Ref. and can 
be expressed in terms of loop integrals, 

[Fl{Xgq), F2{Xgq), F2{Xgq), FiiXgq) , Fi{Xgq)\ 
fc2 



dk 

k'^m? 



(fc2 -I- m^){k'^ -f m^) 



k mj 



2(A:2 + TO?)3' 



2(A:2+to2)3' (P+^2)2' (fc2+^2)5 



(25) 



For later purpose, we give the asymptotic behavior of 

gi{Xgq), 

- Q[9i{x), g2{x), g3{x),g4{x)] 

1, 3 Inx, 3, 6 lux] /x^, for a; 1, 
"1/10, 3/20, 1/2, 1/2], for X ^ 1, (26) 
'1, 1/2, -einx, 1], for a; < 1. 



From Eqs. (^3|-^5|), it is clear that g2A{xgq) correspond 
to photon and gluon attached to the internal squark line, 
while gi,z{xgq) correspond to the opposite case but only 
with the gluon attachment. Note that in the large Nc 
limit, G2(G), 2G2(i?) ^ Nc, while [2G2(i?) - G2(G)] ^ 
0{l/Nc) is suppressed. One always has an Nc factor in 



G^'\ while for C'g\ one only has the Nc factor when a 
gluon attaches to the internal gluino line, which can be 
easily understood by using 't Hooft's double line nota- 
tion. This is also true for the Fi vertex term. How- 
ever, the chiral enhancement factor mg/mb accompany- 
ing 6lr,rl is a unique feature of the F2-term. The mech- 
anism is generic and has been discussed in Ref. |33] , but 
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FIG. 4. Lower and upper limits on squark mixing angles 
Sll,rr obtained by (a) BrSUSY(^+^-)/BrSM(^+^--) > iq%^ 
(b) Br^^^^{p°'y)/Br^^{p°'y) < 4, respectively. 
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FIG. 5. Lower and upper limits on squark mixing angles 
Slr,rl obtained by (a) BrSUSY(^+^-)/BrSM(^+^-) > io%, 
(b) BrSUSY(^o^-)/Bj.SM(.^o^) ^ 4^ respectively. 



SUSY with LR squark mixings gives a beautiful example 

For a direct destructive interference to cut down by 
half the predicted SM rate, one needs the SUSY ampli- 
tude to be 30% of the SM amphtude. This wiU be the 
minimum requirement on SUSY contribution. In Figs. 
|(a) and ||(a) we show limits on \Sll.rr\ and \Slr.,rl\, 
respectively. We require the SUSY contribution alone to 
give 10% of Bi^^{B tttt), corresponding to ~30% in 
amplitude. If we change the required rate contribution 
by a factor k, the values shown in the plots scale by a 
factor y^. From Fig. ^a), we see that the decay rate 
is insensitive to left-left and right-right squark mixings, 
which means gluino box and Sll.rr. related gluino pen- 
guins do not give large contributions. 

The SUSY contribution from gluino box and the -Fi 
term is dominated by Bi^2{xgq), P{xgq) and gi{xgq). For 
m <C mg, B2{xgq) is dominant and contributes through 
cg, which is from gluino box containing dL-bL and qn 
squark lines. For rrig <ti m, P{xgq) is dominant and con- 
tributes through dfi. However, from Eqs. ([l^), (19), 
(p2|), (Q) and (^), we see that this Fi term always re- 
ceive cancellation from Bi{xgq) and gi(xgq), which are 
not too small in this region. Therefore the rise in the 
lower right corner of Fig. ^a) shows insensitivity to 
Sll,rr. as a consequence of this cancellation effect. 

In Fig. ||(a), we show the required \5lr,rl\ to produce 
large enough SUSY contribution in _B — > tttt decay. Note 
that in Eq. (p4[), we use a running mb(/igusY)- For most 



of the parameter space a less than 2% mixing angle in 
left-right mixing is enough to generate such a large SUSY 
contribution. The sensitivity is greatly enhanceded from 
the previous case due to the chiral enhancement factor 
nig/rrif,. Note that there is nothing peculiar about chiral 
enhancement. It only reflects the chiral suppression of 
Cg^j in the SM due to the V— A nature of weak interac- 
tion, which need not be obeyed by interactions beyond 
the SM. 

For left-right mixing case as shown in Fig. ||(a), we 
see that the SUSY contribution is larger for squark mass 
greater than gluino mass and vice versa. For the case of 
heavy squark and light gluino {xgq = m'~/m'^ < 1), the 
gluon preferably radiates off the gluino rather than the 
squark line, as is clear from the behavior of g^{xgq) in 
Eqs. (pJ-Eq). Note that this is the one with Nc enhance- 
ment and therefore gives larger contribution compared to 
TO < rrig case, where the dominant diagrams do not have 
Nc enhancement. The SUSY contribution is dominated 
by 5LR,RLmgg2,{xgq)/{mbm^). For m < m^, as shown 
in Eq. (^), this term becomes SLR,RL'fn^ / C^rn^rnb) and 
is consistent with the sharp rise in the upper left corner 
of Fig. ^(a). The SUSY contribution is small and in- 
sensitive to squark mixing angle in this region. This is a 
generic feature for gluino penguin contributions, as one 
can see from Eq. ( P6[) that all gi/rh? receive a m? /m^ 
suppresion factor in this region. In the reversed case of 
TOg ^ TO, there is only Ijw? suppression, while g^[xgq) 
contribution receive \nxgq enhancement. 

The gluino exchange induced photonic penguin is 
closely related to the gluonic penguin. For example, 
they have similar chiral enhancement behavior as well as 
asymptotic behavior. Recently, the Belle Collaboration 
reports a 90% upper limit on Br(i3° p°7) < 1.06x 10~^ 
, nominally ~ 5 times the SM prediction, but a factor 
of 2 above their previous result reported at ICHEP2000 
1^ . We require the decay rate due to SUSY contribu- 
tion alone to be smaller than 4 times the SM prediction. 



The bounds correspond to 



2|C™|. Note 



that in the LR case, one may have cancellation between 
SM and SUSY . For a direct cancellation, the bounds 
on 5lr^ can be relaxed by 50%. In Fig. ^(b) and ||(b), we 
show the limits on \5ll,rr\ and \5lr.rl\, respectively. 
Similar to the B — > tttt case, the decay rate is insensitive 
to \5ll,rr\, but very sensitive to \5lr,rl\, as expected. 
We also see a sharp rise in the upper left corners of Fig. 
|(b) and|(b), which are similar to ^(a) and indicate that 
SUSY contributions are insensitive to the corresponding 
squark mixings in that region as we have explained in the 
previous paragraph. 

We note that in most of the parameter space shown in 
Fig. |(b), \5lr., alI are constrained to be less than 2%. 
When compared to Fig. |^(a), in most of the parameter 
space \5lr,rl\ impacts more on i? — > p7 than B tttt. 
This can be easily understood by noting that the for- 
mer is a pure loop process while the latter is dominated 
by tree diagrams in SM. Therefore, it is easier for New 
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FIG. 6. Dashed [sohd] lines are upper [lower] bounds 
on squark mixing angles 5lr,rl obtained by Br^u^^(p°7)/ 
BrS^(p''7) < 4 [Br^uSY ^^+^_^^^^sM(^+^--, ^ ^.^^^ 

nig = 200, 500, 700 GeV, respectively. Shaded regions are 
allowed parameter space. 



Physics to affect the former process. 

It is quite interesting that, for the parameter space of 
TO - 300-1000 GeV and m~g < 700 GeV, with mixing an- 
gle ~ 0.2%-0.8%, the model gives sizable contribution to 
B TTTT decay that can account for the smallness of the 
rate, but still satisfy the B ^ p'y constraint. As shown 
in Fig. ^, where the dashed (solid) lines correspond to 
upper (lower) limits on |(5Lfl,flL| from i? — > p7 (tttt), with 
nig = 200, 500, 700 GeV, respectively. The shaded re- 
gions are the allowed parameter space for given nig. For 
nig > 700 GeV we need to > 1 TeV to have allowed 
region, which is beyond the plot. In addition, one can 
also use a smaller (g^), such as to^/4, to enhance the 
color dipole contribution and thus reduce the limit by 
33% and enlarge the overlapping parameter space be- 
tween Figs. ||(a) and ^(b). The existence of this over- 
lap region is closely related to the behavior of gi{xgq). 
From Eqs. (|2^) and (p^), it is clear that for left-right 
mixing, C'^ ccnig/mb giixgq), while C'g (xnig/nib gsixgg). 
For Xgq < 1, g3{x~gq) / gi{xgq) &\\nx~gq\ can be rather 
sizable. Therefore, gluino penguin can give larger contri- 
bution in 6 — > df; than in 6 — *■ ^7 process. 

For illustration, we pick some points from Fig. and 
study the impact of SUSY contributions on Br(7r^7r") 




Jl/2 7C 37t/2 2JI 
arg( 6rl ) 

Br(i3 nn) obtained by using m = 800 GeV, 
and (a) \Slr\ = 0.0035, (b) \Srl\ = 0.0035, 
respectively. The upper band corresponds to the SM predic- 
tion, while the lower band corresponds to the experimental 
result with 2a error range. 



nig = 200 GeV, 





Jl/2 71 3JI/ 2 271 71/2 n 3n/2 2ji 

arg(5LR,RL) arg( 5lr, rl ) 

FIG. 8. (a) Upper (lower) lines corespond to Br(p''7)/ 
Br^'^(p"7) with rrig =200 (500) GeV, m=800 (900) GeV, 
15LB,flL 1=0.0035 (0.006). (b) The asymmetry in B -» tt+tt" 
with same parameter space as case (a). 



and Br(p"7). In Fig. we show Br(7r+7r ) obtained by 
using m~g = 200 GeV, to = 800 GeV, and (a) \5lr\ = 
0.0035, Srl = Sll,rr = 0, (b) \Srl\ = 0.0035, Slr = 
Sll,rr. — 0, respectively. The upper band corresponds 
to the SM prediction, while the lower band corresponds 
to the averaged experimental result Br(7r+7r^) — 4.4 ± 
0.9 with 2(7 error range. With a,rg{6LR,RL) within the 
dashed hnes, i.e. a.ig{6LR) ~ 4.3-27r, a,rg{6RL) ~ 1.2-3.2, 
Br(7r+7r^) can be brought down by SUSY contributions 
to experimental range. The strength factor of ap^ ^3 36 



sin 29 = 



2\C^C'\ 



(27) 



can be as large as 90% in this case. The measurability 
of the asymmetry in i? — > p7 decay is better than in 
B ^ K*j since it readily provides vertex information 
|6|. 

It is clear from Fig. |^ that we might as well take 
nig = 500 GeV, to = 900 and Slr,rl = 0.006. The pre- 
vious case corresponds to Xgq — 0.06, while in this case 
we have larger Xgq — 0.31. These two cases also represent 
other cases with similar Xgq, while to need not be that 
heavy. The figure for Br (tt+tt") are almost identical to 
Fig. 1^. However, as we show in Fig. ^(a), the latter case 
has greater Bi{p'^^). Note that RL case is insensitive to 
arg(^i?L), since the rate is proportional to |C-yp -I- \C!y\'^. 
In this case, Br(p"7)= Br^'^(/9°7)-FBrSUSY|-^o^) jg within 
5 times of SM rate as required by Fig. ^ and Fig. ^. All 
B ^ TTTT favored range, arg{6RL) ^ 1.2 - 3.2 are also al- 
lowed by the B ^ pj constraint. The sin 26* is 80%. For 
the LR case, the induced C^^^^ may have constructive or 
destructive interference with as arg(<5L_R) changes. 

Within the a.rg{5LR) ~ 4.3-2tt range, allowed by i? ^ tttt 
rate, there are quite some parameter space to satisfy the 
B ^ pj constraint. The rate can be close to the SM ex- 
pectation. In Fig. H(b), we show the asymmetry, a^+^-. 
Note that the SUSY prediction for ajr+Tr- from the pre- 
vious two paramter points are close as in the Br(TT"''Tr~) 
case. In generalized factorization, _ ~ 3.3% with 
our input parameter, a smaller {q ) will have a slightly 
larger asymmetry. With SUSY, a^+T^- can be rang- 
ing within —1% to 8%. Note that in QCD factoriza- 
tion approach including weak annihilation contribution 



-5% to 15%, for 



60° II6 . It is difhcuh to 



8 



distingiush SUSY contributions from the SM prediction 
from a_+^- . 



IV. DISCUSSION 

Left-left and/or right-right Q d-b mixings with few % 
to few 10% mixing angle can generate large enough con- 
tribution to _Bc;-mixing that could reduce aj/^j^^ from 
its SM value. As shown in Fig. ^b), such squark mixing 
angles are safe from _B — > p7 constraint. However, as 
shown in Fig. ^(a), one needs large mixing with sizable 
mass splitting to affect i? — > tttt decay rate in this case. 
Such a large mixing is already ruled out by the exper- 
imental measurement of Atob^, as shown in Figs. |l|(a) 
and |^(b), unless one fine tunes the parameter space to 
be very close to Xgq = 2.43, and turn off left-left or right- 
right mixings. In other words, one needs high degree of 
fine tuning to account for both low sin 20i and B — + tttt 
decay rate with LL or RR mixings. It is much easier to 
compete with the SM box diagram and modify sin 2(j)i 
than to compete with tree dominated B —^ tttt decay. 

Alternatively, left-right and/or right-left d-b mixings 
with few % mixing angles could also give sizable con- 
tribution to Bd mixing. However, because of the am- 
plification effect of chiral enhancement, the size of this 
mixing angle is severely constrained hy B —f pj, to be 
less than 2% in most of the parameter space given in 
Fig. ||(b). It cannot be the source that gives sizable 
contribution to B^ mixing, as one can tell by comparing 
Figs. H and||(b). It is interesting that, as noted already 
in the previous section, there is parameter space where 
rrig is suitably light and the mixing angle 6lr,rl is less 
than 1%, where the model gives sizable contribution to 
B —^ TTTT decay without violating B pj constraint. In 
other words, we need left-right and/or right-left mixings 
rather than left-left and/or right-right mixings to affect 
B ^ TTTT decay. Thus, if the smallness of Br(i3 tttt) is 
due to SUSY, it is likely that one will have large effects in 
b ^ d'j, including rate enhancement and mixing induced 
asymmetry |3^, which can be easily close to 100%. 

In Section III, we used 4>3 = 65°, which is a CKM fit 
like value, since SUSY contribution to _B — > tttt decay is 
uncorrelated with the SUSY contribution to B^ mixing. 
The CKM fit may still be viable and need not support 
large (ps as a solution of low B ^ tttt rate. 

It is clear that we still need correct interference pat- 
terns, i.e. correct SUSY phases, to reduce aj/^p^^ and 
Br(i3" 7r+7r~). Since the effects arise from differ- 
ent squark mixing sources, one can always find separate 
SUSY phases to achieve this. As we show in Section II 
and Section HI, we may have accessible allowed region 
on SUSY phases. 



V. CONCLUSIONS 

We have shown that it is possible for SUSY models to 
account for the smaller sin2(/)i and Br(i? — > tttt) values 
that seem to be emerging from the B Factories. However, 
they would have to come from different flavor mixing 
sources. The smallness of sin 201 is most likely arising 
from left-left (right-right) squark mixing, while the deficit 
in Br(i? tttt) is most likely due to left-right squark 
mixings. The two are basically uncorrelated. 

Because of similarity in chiral enhancement, the loop 
induced b ^ dj process is even more sensitive to Slr.rl 
than the tree dominated B —f tttt decay. Therefore, if 
SUSY affects the latter, the effects in the former would 
be even more prominent. We emphasize that _B — > p7 
could be considerably larger than expected in SM if the 
smallness of B ^ tttt rate is in part due to SUSY. 
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